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ABSTRACT 

There are numerous forkhead transcription factors in 
mammalian cells but we know little about the molec- 
ular functions of the majority of these. F0XK2 is 
a ubiquitously expressed family member suggest- 
ing an important function across multiple cell types. 
Here, we show that F0XK2 binds to the SIN3A and 
PR-DUB complexes. The PR-DUB complex contains 
the important tumour suppressor protein, the deu- 
biquitinase BAP1. F0XK2 recruits BAP1 to DNA, 
promotes local histone deubiquitination and causes 
changes in target gene activity. Our results therefore 
provide an important link between BAP1 and the tran- 
scription factor F0XK2 and demonstrate how BAP1 
can be recruited to specific regulatory loci. 

INTRODUCTION 

The forkhead transcription factor family contains over 40 
members and all of these contain a forkhead winged helix- 
turn-helix DNA binding domain (1,2). This domain is re- 
sponsible for sequence-specific DNA binding, generally to 
sites resembling the RYMAAYA (R = A or G; Y = C or 
T; M = A or C) core motif ((3,4); reviewed in (1)). Many 
of these proteins are expressed in a restricted or cell type- 
specific manner but a group of these are ubiquitously ex- 
pressed across many cell types, including the related pro- 
teins FOXKl and FOXK2. In addition to the forkhead 
DNA binding domain, FOXKl and FOXK2 also contain 
a FHA domain in their N-terminal regions. FOXKl has 
been the best studied of these proteins and the mouse ho- 
mologue of FOXKl, MNF, has been shown to regulate 
the proliferation of myogenic stem cells (5,6). Molecularly, 
FOXKl /MNF acts in a repressive manner and can do so 



through recruitment of SIN3 -containing co-repressor com- 
plexes (7,8). Human FOXKl has been less studied but can 
associate with serum response factor (SRF) and modulate 
its transcriptional activity (9). Comparatively less is known 
about the function of FOXK2. FOXK2 was initially iden- 
tified as a regulator of IL-2 transcription (10). Molecularly, 
few binding partners have been identified but it can bind to 
transforming proteins adenoviral El A and papillomavirus 
E6 (11) and to API family members (12). Little is known 
about upstream regulatory pathways but a hnk to the cell 
cycle is suggested by the observation that FOXK2 is phos- 
phorylated by cyclin dependent kinase (CDK)-cychn com- 
plexes (13). In addition to potential transcriptional roles, an 
alternative function in DNA repair is suggested by the ob- 
servation that FOXK2 can bind to DNA containing G/T- 
mismatches (14). 

Our genome-wide ChlP-seq studies showed that FOXK2 
binds to thousands of regions spread throughout the 
genome where it can promote API binding at a subset of 
sites (12). One possible role for FOXK2 in this context 
would be to modify chromatin structure but how FOXK2 
might achieve this is unclear. Here, we have searched for co- 
regulatory partners for FOXK2 and identified the SIN3 and 
polycomb repressive deubiquitinase (PR-DUB) complexes. 
We focussed on verifying functional interactions with the 
PR-DUB complex as this contains the important tumour 
suppressor protein BAPl (reviewed in (15)) and demon- 
strated that FOXK2 acts as a sequence- specific binding 
transcription factor to recruit BAPl to chromatin and pro- 
mote local histone deubiquitination. 

MATERIALS AND METHODS 

Plasmid constructs 

The following plasmids were used in mammalian cell trans- 
fections. pAS 1178 (pCDNA3-GAL4-FOXK2-FH A), 



To whom correspondence should be addressed. Tel: +44 161 275 5979; Fax: +44 161 275 5082; Email: a.d.sharrocks@manchester.ac.uk 
Present address: 

Namshik Han, Gurdon Institute and Department of Pathology, University of Cambridge, Tennis Court Road, Cambridge CB2 IQN, UK. 
© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.Org/licenses/by/3.0/), which 
permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. 



Nucleic Acids Research, 2014, Vol 42, No. 10 6233 



has been described previously. The following plasmids 
were constructed to express N-terminally Flag-tagged 
human BAPl: pAS1443 [pCDNA3-Flag-BAP 1(1-729)], 
pAS1447 [pCDNA3-Flag-BAPl(l-598)] and pAS1448 
[pCDNA3-Flag-BAPl(241-729)], which were generated 
by ligating BamHI/NotI double digested PCR products 
generated with the oHgonucleotide pairs ADS3817/3819, 
ADS3817/3821 and ADS3820/3819, respectively, on the 
template pOTB7-BAPl (GenomeCube, IMAGE ID: 
3543914), into the same sites from the vector pAS3010. 
pAS3010 contains F0X03 cloned in frame with an N- 
terminal Flag tag (using a BamHI site) into pCDNA3 
(kindly provided by Alan Whitmarsh). 

The plasmids used in glutathione ^-transferase 
(GST) pulldown experiments were pAS2194 [encoding 
GST-FOXK2(l-218) and pAS2195 (encoding GST- 
FOXK2( 189-660))] (13). The FOXK2 FHA domain mu- 
tant GST-FOXK2(l-218)(R58A) (pAS1450), was made 
by Quikchange mutagenesis using the primer-template 
combination ADS3825/ADS3826-pAS2194. 

Tissue culture, cell transfection and RNA interference 

U20S and HEK293T cells were grown in Dulbecco's mod- 
ified Eagle's medium (DMEM) supplemented with 10% 
foetal bovine serum. U20S cells stably expressing FOXK2- 
HF (U20S-FOXK2-HF) (12) were grown in DMEM 
supplemented with 1.5 |xg/ml puromycin and 10% foetal 
bovine serum. 

To construct U20S cell lines stably expressing shFOXK2 
or a matched control cell line, non-silencing pGIPZ lentivi- 
ral shRNAmir control vector (RHS4346) and human 
GIPZ FOXK2 shRNA (RHS4430-101517573, clone ID: 
V3LHS_308724) were transfected into U20S cells using 
TurboFect Transfection Reagent (Thermo Scientific). Cells 
were then subjected to selective pressure using 5 |xg/ml 
puromycin. Single cells were then expanded into colonies 
to generate individual clonal cell lines. FOXK2 protein level 
was barely detectable by immunoblotting in clonal lines 5-1 
and F5-1 that were subsequently used in this study. 

For plasmids transfection, 3x10^ HEK293T cells were 
transfected with 8 [xg of plasmid. Transfections were carried 
out using Polyfect (Qiagen) transfection reagent. After 24 h, 
cells were harvested and used for co-immunoprecipitation 
assays. 

siRNAs for FOXK2 and BAPl were ON-TARGETplus 
SMART pools from Dharmacon (L-008354-00-0005 and 
L-00579 1-00-0005), control non-targeting siRNAs (Dhar- 
macon, D-00 18 10-0 1-50) were used throughout. To carry 
out RNA interference (RNAi), cells were transfected with 
50 nM siRNA using Lipofectamine® RNAiMAX (Invit- 
rogen) according to the manufacturer's instructions. Forty- 
eight hours after transfection, the cells were harvested for 
further analyses. 

Western blot analysis, co-immunoprecipitation and GST 
pulldown assays 

Western blotting was carried out with the primary anti- 
bodies; FOXK2 (A301-729A), HCFCl (A301-400A) and 
ASXL2 (A302-037A) from Bethyl Laboratories; BAPl (C- 
4, SC-28383), mSIN3A (K20, sc-994), Gal4 (sc-577) and 



ERK2 (sc-154) from Santa Cruz Biotechnology; Flag M2 
(F3165) from Sigma; HDACl(abl767) from Abeam. West- 
ern blots derived from whole cell lysates or immunoprecipi- 
tated proteins were visualized after incubation with primary 
antibodies using Infrared IRDye-labelled secondary anti- 
bodies and the signal was collected with a Li-Cor Odyssey 
infrared imager. 

Co-immunoprecipitation was performed as previously 
described (12), antibodies used for immunoprecipitation: 
FOXK2 (anti-ILFl antibody, ab5298. Abeam); antibodies 
for Flag, BAPl and Gal4 as described above; for each im- 
munoprecipitation, normal IgG was used as negative con- 
trol: rabbit (12-370), mouse (12-371) IgG from Millipore 
and normal goat IgG (sc-2028) from Santa Cruz Biotech- 
nology. For interrogating phosphorylation-dependent in- 
teractions, co-immunoprecipitated proteins were treated 
with 400 units of X-phosphatase (NEB P0753S) for 30 min 
at 30°C followed by four further washes in immunoprecipi- 
tation buffer prior to loading on a gel. 

GST pulldown analysis was performed essentially as de- 
scribed previously (16), with purified bacterially expressed 
recombinant GST-FOXK2 derivatives and either total cell 
extracts from U20S cells or in vitro translated ^^S-labelled 
BAPl protein synthesized using the TNT T7 Quick-coupled 
transcription/ translation system (Promega). Where indi- 
cated, ethidium bromide (200 jxg/ml) was added to the GST 
pulldown reactions. 

Proteomic analysis 

Rapid immunoprecipitation mass spectrometry of endoge- 
nous proteins (RIME) was performed as described previ- 
ously (17). U20S-FOXK2-HF cells and the control 'empty 
vector' U20S-HF cells were grown in complete DMEM in 
the presence of 1.5 |JLg/ml puromycin. Anti-Flag M2 anti- 
body (F3165, Sigma) was used for immunoprecipitation of 
cross-linked FOXK2. 

The interaction networks of FOXK2 associated proteins 
identified using RIME were analysed using STRING (18), 
with parameters set for Active Prediction Methods as ex- 
periments, textmining and databases. 

Immunofluorescence analysis and proximity ligation assay 

Immunofluorescence experiments were performed as de- 
scribed previously (13). DNA was stained with 4',6- 
diamidino-2-phenylindole (DAPI) using mounting medium 
of ProLong Gold Antifade Reagent (Invitrogen). The 
Duolink In Situ Red Starter Kit Mouse/Rabbit (Sigma) was 
used for proximity hgation assay (PL A) experiments. U20,S 
U20S-FOXK2-HF or shFOXK2 clonal line 5-1 cells were 
seeded onto coverslips at 5 x 10^ cells per well in a 12- 
well plate. The cells were fixed with 4% paraformaldehyde- 
phosphate buffered saline (PBS) at room temperature for 
30 min and permeabilized using 0.2% Triton-XlOO-PBS 
for 6 min. After permeabilization, the cells were incubated 
in the blocking buffer (5% bovine serum albumin (BSA)- 
PBS, freshly prepared) for 30 min at room temperature 
(RT). Cells were then incubated with the primary antibod- 
ies diluted in the antibody diluents for 1 h at RT. For the 
rest of the assay, the manufacturer's instructions were fol- 
lowed. The coverslips were mounted using DUOLINK in 
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situ mounting medium with DAPI. Images were collected 
on an Olympus BX51 upright microscope using a x63 oil 
immersion objective and captured using a Coolsnap ES 
camera (Photometries) through MetaVue Software (Molec- 
ular Devices). Images were then processed and analysed us- 
ing Image J. To quantify the FLA signals, 5-10 image sec- 
tions were analysed for the ratio of signals per nuclei. The 
negative control level was set to 1 . 

RT-PCR and chromatin immunoprecipitation (ChIP) assays 

Total RNA was isolated using a RNeasy kit (Qiagen) and 
transcripts detected in a one-step RT-PCR reaction using 
Quantitect SYBR green reagent (Qiagen). The primer-pairs 
used for RT-PCR experiments are listed in Supplementary 
Table SI A. 

Chip assays using control IgG (Millipore) or antibodies 
specific to BAPl, FOXK2, HCFCl, ASXL2, Flag, H2A 
(ab 18255, Abeam) and ubiquityl-histone H2A (Lysll9) 
(H2Aub, D27C4, 8240S, from Cell Signaling Technology), 
were performed essentially as described previously (12), us- 
ing 1 X 10^ to 5 X 10^ U20S cells for a standard ChIP 
Bound regions were detected by quantitative PCR (qPCR) 
(using primers listed in Supplementary Table SIB), from at 
least two independent experiments, using Quantitect SYBR 
green PCR reagent (Qiagen). Results were analysed with 
Rotorgene Q software (Qiagen) relative to input using the 
standard curve method. ChlP-seq analysis for Flag-tagged 
FOXK2 and Flag-tagged BAPl from U20S-FOXK2-HF 
and U20S-BAP1-HF cells was performed essentially as de- 
scribed previously (12). 

Statistical analysis for real-time PCR results was per- 
formed using the Student's ^-test. The error bars in all 
graphs represent standard deviation. 

Bioinformatics analysis 

For microarray analysis, we first examined our previously 
published datasets following treatment of U20S cells with 
either siFOXK2 or siGAPDH as a control (Array Express 
E-MEXP-3106) (12). Each experiment was performed us- 
ing the Affymetrix hgul33plus2.0 microarray platform, and 
repeated in triplicate. Raw intensity CEL files for these 
six datasets were normalized, and converted to estimates 
of probe set expression level, using the robust multiar- 
ray averaging (RMA) algorithm from the affy (19) pack- 
age in R. Differential expression between siFOXK2 and 
siGAPDH treated cells was then estimated using unpaired 
tests within the limma package in R (20). Probe set iden- 
tifiers were converted to protein-coding gene symbols (En- 
sembl V.73 gene symbols, and HGNC identifiers) using the 
most recent mappings provided by Affymetrix. In all con- 
texts, we used the Ensembl v.73 protein-coding gene collec- 
tion, obtained through the BioMart portal (vO.7), retaining 
both known and novel protein-coding transcripts (21). In 
total, 18 409 Ensembl gene identifiers matched at least one 
probe set on the hgul33plus2.0 platform. We also examined 
externally published datasets on the same microarray plat- 
form (hgul33plus2.0), following treatment of U20S cells 
with either shBAPl or a control shRNA (GEO: GSE23035) 
(22). Two different shRNAs specific to BAPl were used. 



and all experiments repeated in triplicate, giving nine ar- 
rays in total. Raw intensity CEL files were normalized as 
before using RMA, and differential expression between the 
six shBAPl arrays and three shControl arrays estimated us- 
ing the limma package in R. 

For ChlP-seq analysis, we used our Illumina ChlP-seq 
read hbraries for FOXK2-bound compared to input DNA 
in U20S cells (ArrayExpress: E-MTAB-2204), together 
with externally published ChlP-seq read libraries for DNA 
bound by HCFCl, compared to input, in HeLa cells (GEO: 
GSE31417) (23). We did not reanalyse previously pubhshed 
SOLID ChlP-seq read libraries (12) due to more recent find- 
ings that this technology suffers from bias against CpG- 
enriched regions, including many active gene promoter re- 
gions (our unpublished data). Read libraries were mapped 
to the human genome (hgl9) using Bowtie vl.O with strin- 
gent settings (v = 1, m = 0) (24). In all cases, peak-calling 
was carried out using HOMER (v4.2), searching for narrow 
binding regions (style = factor), and with default settings 
(Poisson threshold p = 10~^) (25). Peak locations were com- 
pared to one another and to protein-coding gene transcrip- 
tion start sites (TSSs) using in house methods written in Java 
(code available upon request). The density of FOXK2 and 
HCFCl tags within 5 kb regions centred on FOXK2 peak 
centres was visualized using Seqminer vl.3.3 (26). 

Statistical analysis 

To annotate genes as up-regulated, unchanged, or down- 
regulated, upon siFOXK2 or siBAPl treatment, we im- 
posed a range of fold change (fc) cut-offs (1.2, 1.3, 1.4 and 
1 .5) compared to control siRNA treatment (^nominal < 0.05). 
Rarely, a gene with multiple probe sets passed thresholds 
for both up- and down-regulation, and in such cases the 
genes were omitted from further statistical analysis. At each 
fold change cut-off, we then measured the significance of 
overlaps between gene lists up- or down-regulated upon ei- 
ther siFOXK2 or shBAPl using the hypergeometric test. We 
found results were independent of the fold change param- 
eter (data not shown), so we report the results for the rel- 
atively less stringent /c = 1.2 (up-regulation) or /c < 1/1.2 
(down-regulation) . 

RESULTS 

FOXK2 binds to the BAPl-containing PR-DUB complex 

To further our understanding of FOXK2 function, we used 
RIME (17) to identify novel FOXK2 binding partners. This 
technique is designed to maximize the chances of isolating 
factors that are associated with chromatin-bound FOXK2, 
by using formaldehyde-mediated cross-linking prior to im- 
munoprecipitating FOXK2 (Figure lA). Immunoprecipta- 
tions (IPs) were performed using anti-Flag antibody using 
extracts from cross-linked U20S-FOXK2-HF cells which 
harbour Flag-tagged FOXK2 expressed at endogenous lev- 
els (12), and the resulting immunoprecipitates were anal- 
ysed by mass spectrometry. This resulted in the identifica- 
tion of FOXK2 as expected (32 peptides giving 65% cov- 
erage) but also 119 proteins that were not also identified 
in control IPs from U20S-HF cells. Amongst these pro- 
teins were proteins belonging to the SIN3A core complex 
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Figure 1. Identification of FOXK2 binding proteins by RIME. (A) Schematic illustration of the RIME protocol used to identify FOXK2 associated factors 
in the context of chromatin association. (B) Summary of FOXK2 interaction proteins belong to the SIN3A core complex and the PR-DUB complex. 
HCFCl is shown in brackets as it is unclear whether this is part of the core PR-DUB complex. (C) Visualization using STRING of a sub-network of 
interactions between FOXK2 binding proteins. (D) Validation of FOXK2 interaction with endogenous components of SIN3 A and PR-DUB core complex 
proteins by co-immunoprecipitation (IP) experiments using anti-Flag (FOXK2) antibody in U20S-FOXK2-HF cells. Precipitated proteins were detected 
by immunoblotting (IB) using the antibodies as indicated. Arrows represent the bands corresponding to each of the full-length proteins. Note that the 
ASXL2 blot is from a different IP experiment. (E) PLA analysis of interaction between Flag-tagged FOXK2 and endogenous HCFCl in U20S-FOXK2- 
HF (top and bottom left panels) or U20S (bottom right) cells. The combinations of antibodies used are shown above and below each panel (IgG represents 
a non-specific antibody). DNA is stained using DAPI (blue) and the PLA signal is red. Quantitative analysis of PLA signals in the nucleus is shown below. 
The level of signals/nuclei in the control PLA sample (Flag and non-specific IgG antibodies) was set as 1. 



and the PR-DUB complex (Figure IB). The association 
with a SIN3-hke complex was expected as the closely re- 
lated protein FOXKl has previously been shown to bind 
to SIN3A and SIN3B (7,8) whereas binding of FOXK2 to 
the PR-DUB complex was novel. Interestingly, several of 
the proteins in these two complexes were also identified in 
reciprocal IPs for BAPl interacting proteins (Figure IB) 
(27,28). Moreover, analysis of the FOXK2 interactome us- 
ing STRING (18) revealed an interconnected sub-network 
containing both the SIN3A and PR-DUB complexes, with 
binding to HCFCl providing the links between the two 
complexes (Figure IC). Interactions between Flag-tagged 
FOXK2 and endogenous HDACl and SIN3A components 
of the SIN3A core complex were verified by co-IP (Fig- 
ure ID). Similarly, we were able to verify FOXK2 interac- 
tions with the PR-DUB complex components BAPl and 
ASXLl but were only able to detect weaker co-associations 
with HCFCl in co-IP experiments (Figure ID). To establish 
whether HCFCl interacted with FOXK2 in vivo, we instead 
used the in situ PLA which can detect protein-protein in- 
teractions in cells with high specificity and sensitivity. Incu- 



bation of cells stably expressing Flag-tagged FOXK2 with 
either the Flag or HCFCl antibody alone resulted in only 
background signals. However, co-incubation with both an- 
tibodies together gave a strong PLA signal but this was only 
apparent when the cells contained Flag-tagged FOXK2 and 
was not observed in parental U20S cells. These results are 
indicative of a close association of HCFCl and FOXK2 
(Figure IE). 

BAPl has recently been implicated as an important tu- 
mour suppressor protein in a number of different tumour 
types (reviewed in (29)) and found to be one of seven genes 
whose mutation is the most indicative of poor prognos- 
tic outcome (30). We therefore focussed on the PR-DUB 
complex and on BAPl in particular for further investiga- 
tion. First, we demonstrated that endogenous FOXK2 and 
BAPl interact in reciprocal co-IP experiments (Figure 2 A 
and B). Both proteins are found in the nuclear pellet frac- 
tion (Supplementary Figure SI) and immunofluorescence 
microscopy showed that their expression patterns within 
the nucleus show substantial overlap (Figure 2C). To gain 
a higher resolution view of their co-localization, we again 
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Figure 2. Interactions between endogenous FOXK2 and BAPl. Re- 
ciprocal co-immunoprecipitation experiments using either FOXK2 (A) 
or BAPl (B) antibodies for immunoprecipitation (IP) from U20S cells. 
Co-precipitated endogenous BAPl and FOXK2 were detected by im- 
munoblotting (IB). Arrows represent the bands corresponding to each of 
the full-length proteins. (C) Immunofluorescence detection of endogenous 
FOXK2 (green) and BAPl (red) and their co-localization in the nucleus 
[indicated by DAPI staining of DNA (blue)]. (D) Imaging and quantifica- 
tion of PLA signals generated by the indicated combinations of antibod- 
ies (IgG represents a non-specific antibody). DNA is stained using DAPI 
(blue) and the PLA signal is red. Quantitative analysis of PLA signals in 
the nucleus is shown on the right. The level of signals/nuclei in the control 
PLA sample (BAPl and non-specific IgG antibodies) was set as 1. 



turned to the PLA assay to examine whether FOXK2 and 
BAPl are closely juxtaposed in vivo. Incubation of cells 
with either the FOXK2 or BAPl antibody alone resulted 
in only background signals. However, co-incubation with 
both antibodies together gave a strong PLA signal (Figure 
2D, left panels). Importantly, this signal was equivalent to 
that seen between BAPl and either the PR-DUB compo- 
nent ASXL2 or the closely associated protein HCFCl (Fig- 
ure 2D, right panels). Importantly the PLA signals between 
FOXK2 and BAPl were lost in a stable cell line where en- 
dogenous FOXK2 is depleted by expression of an shRNA 
against FOXK2 (Figure 2D, bottom right panel). Together 
these results firmly establish BAPl from the PR-DUB com- 
plex as an interaction partner for FOXK2. 



Mapping the FOXK2-BAP1 interaction domains 

To provide insights into complex formation between 
FOXK2 and BAPl, we first mapped the region of FOXK2 
responsible for binding to BAPl. First we carried out in 
vitro GST pulldown assays with GST fusions to the N- 
and C-terminal parts of FOXK2 (Figure 3A). Strong in- 
teractions with BAPl were seen with FOXK2(l-218) con- 
taining the N-terminal region whereas background bind- 
ing was seen to the C-terminal FOXK2( 189-660) construct 
(Figure 3B, lanes 4 and 5). Binding of BAPl to the N- 
terminal region of FOXK2 in vivo was confirmed by co- 
IP of a GAL-FOXK2(l-218) fusion protein (Figure 3C). 
To rule out the possibihty that the interactions we de- 
tect are indirect and mediated by contaminating DNA, we 
performed a GST pulldown assay with GST-FOXK2(l- 
218) and total U20S cell ly sates and investigated BAPl 
interactions in the presence and absence of ethidium bro- 
mide. The inclusion of ethidium bromide had no effect 
on FOXK2-BAP1 interactions (Figure 3D), demonstrating 
that the interactions are not DNA-mediated. FOXK2(l- 
218) contains the FHA domain which has been shown to 
be a phosphopeptide binding domain in other proteins (re- 
viewed in (31)). Further truncations of the regions located 
N- and C-terminally to the FHA domain created proteins 
FOXK2(52-218) and FOXK2(l-128) which were still ca- 
pable of binding FOXK2 in vitro (Supplementary Figure 
S2) implicating the FHA domain as important for BAPl 
binding. However, the FHA domain alone was insufficient 
for BAPl binding, possibly due to incorrect folding in this 
truncated construct. Given the fact that the FHA domain 
is important for BAPl binding, we tested an R58A mu- 
tant version of FOXK2 which is predicted to ablate its 
phosphopeptide binding activity (31). Reduced binding of 
BAPl was observed with FOXK2(R58A) (Figure 3E, top 
panel). Moreover, binding to both HCFCl and SIN3A was 
also reduced. This suggested a role for phosphorylation- 
dependent interactions between FOXK2 and its interaction 
partners. To test whether this was indeed the case, we re- 
peated a co-immunoprecipitation experiment, but treated 
the co-precipitates with X phosphatase to remove any phos- 
phate groups. Both BAPl and FOXK2 showed mobility 
changes, indicative of dephosphorylation, but there was no 
effect on FOXK2-BAP1 interactions (Figure 3F). Further- 
more, treatment of cells with the CDK inhibitor alster- 
paullone had little effect on FOXK2-BAP1 interactions 
(Supplementary Figure S2C). It therefore appears likely 
that the R58A mutation does more than affect phospho- 
dependent interactions but underscores the importance of 
this domain in mediating interactions between FOXK2 and 
BAPl. 

Next, we turned to BAPl and deleted either the N- 
terminal ubiquitin carboxyl-terminal hydrolase (UCH) do- 
main or the C-terminal domain (CTD) which contains the 
UCH 3 7 -like domain (ULD) and nuclear localization signal 
(NLS) (Figure 3G). Surprisingly, deletion of either domain 
reduced interactions between Gal-FOXK2(l-218) in co-IP 
experiments, although the lower expression of BAP 1(241 - 
729) likely contributed to the lower levels of interaction 
observed (Figure 3H, first and third panels). Similarly, in- 
teractions between BAPl and HCFCl were reduced upon 
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Figure 3. Mapping the FOXK2-BAP1 interaction domains. (A) Schematic representation of full-length FOXK2 and the GST fusions to the C-terminal 
(amino acids 1-218) and N-terminal (amino acids 189-660) regions of FOXK2. The positions of the forkhead associated (FHA) and forkhead domains 
(FOX) are shaded in grey. (B) GST pulldown assays using bacterially expressed GST-tagged FOXK2(l-218) and FOXK2( 189-660) and in vitro translated 
BAPl. Arrows mark the positions of full-length GST-FOXK2 fusion proteins. A total of 10% input is shown. (C) Co-immunoprecipitation experiments 
to analyse interactions between FOXK2(l-218) and BAPl. HEK293T cells were transfected with the indicated plasmids encoding FOXK2(l-218) fused 
with the Gal4 DNA binding domain and Flag-tagged BAPl, followed by immunoprecipitation (IP) with anti-Flag antibody and immunoblotting (IB) 
with either an anti-Flag or an anti-Gal4 antibody. ERK2 is a loading control. The asterisk marks a non-specific signal. (D and E) GST pulldown assay 
using either GST, or wild-type (WT) or FHA mutant (R58A) versions of the GST-FOXK2(l-218) fusion protein and total cell extracts from U20S cells. 
Interacting BAPl, HCFCl, SIN3A and input GST fusion proteins were revealed by IB. Arrows mark the positions of full-length GST-FOXK2 fusion 
proteins. A total of 3% cell lysate input is shown. Ethidium bromide was added to the GST pulldown reactions where indicated. (F) Co-immunoprecipitation 
experiments using FOXK2 antibodies for IP from U20S cells. Co-precipitated endogenous BAPl was detected by IB. Where indicated, the final co-IP was 
treated with X phosphatase. (G) Schematic representation of full-length BAPl and the indicated N- and C-terminal truncation mutants. The locations of 
the ubiquitin carboxyl-terminal hydrolase (UCH) domain and the C-terminal domain (CTD) are shown. (H) Co-immunoprecipitation analysis of FOXK2 
interactions with BAPl deletion mutants. HEK293T cells were transfected with the indicated plasmids encoding FOXK2(l-218) fused to the Gal4 DNA 
binding domain and Flag-tagged full-length or truncated mutants of BAPl, followed by IP with anti-Flag antibody and IB with the anti-Gal4, anti-HCFCl 
and anti-Flag antibodies. The asterisk marks a non-specific signal. 



truncation of either its N- or C-terminal domains (Fig- 
ure 3H, second panel). As loss of the UCH domain in 
the BAPl(241-729) construct reduces FOXK2 binding we 
tested whether the deubiquitination activity of BAPl is re- 
quired for FOXK2 binding. However, the catalytically in- 
active BAP1(C91S) mutant bound with equal efficiency as 
wild-type (WT) BAPl to FOXK2 (Supplementary Figure 
S3). 

Together these results suggest that BAPl has to be struc- 
turally intact to bind to FOXK2 and that binding occurs 
through the N-terminal FHA domain of FOXK2. 



FOXK2 and the PR-DUB bind to the same genomic regions 

Having established that FOXK2 and BAPl interact in vivo, 
we next established whether they could bind to the same 
genomic regions. Previously we identified the FOXK2 ge- 
nomic binding regions in U20S cells by ChlP-seq (12). We 
therefore tested four of these for endogenous BAPl binding 
(Figure 4A) and verified FOXK2 occupancy (Supplemen- 
tary Figure S4) by qPCR-ChIP Binding of both FOXK2 
and BAPl was detected at all of these regions while the con- 
trol intronic region from the MCM3 locus exhibited bind- 
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Figure 4. FOXK2 and the PR-DUB bind to the same chromatin regions. ChlP-qPCR analysis of endogenous BAPl (A) and the other PR-DUB compo- 
nents HCFCl and ASXL2 (B) binding to the indicated FOXK2 genomic binding regions. Blacks bars in (B) represent binding to FOXK2 binding regions 
associated with the MCM3 and KDM3A loci; white bars are negative control region which does not bind to FOXK2 (MCM3-mt9). Data are normalized 
against input DNA and shown relative to binding to non-specific IgG (taken as 1). Data are the average of three independent experiments. (C) Genome 
wide correlation between FOXK2 and HCFCl (23) binding regions. Heatmap of HCFCl and FOXK2 read densities mapped onto FOXK2 peak summits 
and ranked according to FOXK2 signal. (D) UCSC genome browser view of FOXK2 and HCFCl binding profiles associated with the DDX19A and 
VPS51 loci. 



ing levels close to background binding observed with con- 
trol IgG. Binding of Flag-tagged BAPl was also detected at 
several of these regions (Supplementary Figure S5A). Fur- 
thermore, we could also detect binding of the other PR- 
DUB complex member ASXL2 and HCFCl at the FOXK2 
binding regions associated with the MCM3 and KDM3A 
loci (Figure 4B). We also compared the genome- wide bind- 
ing profile for FOXK2 (12) with that of HCFCl (23), and 
found a strong co-association (Figure 4C), despite the dif- 
ferent cell types being studied (U20S versus HeLa). Indeed 
60% (2203/3674) of the summits of the HCFCl binding re- 
gions are located within 200 nucleotides of a FOXK2 bind- 
ing region summit as exemplified by the regions associated 
with the DDX19A and VPS51 loci (Figure 4D). Similar 
ChlPseq datasets are not available for BAPl, and our at- 
tempts to generate such datasets using stable cell lines ex- 



pressing Flag- tagged BAPl yielded only a handful of bind- 
ing regions, therefore we were unable to do a global com- 
parison for BAPl. 

Functional interactions between FOXK2 and BAPl in target 
locus regulation 

To establish whether FOXK2 and BAPl functionally inter- 
act, we first compared two microarray studies performed 
in U20S cells where either FOXK2 (12) or BAPl (22) had 
been depleted by siRNA and shRNA treatment respectively 
to identify commonly regulated genes. A large number of 
genes were identified which were either commonly down- 
regulated (249 genes) or up-regulated (371 genes) upon de- 
pletion of either factor (Figure 5A). Importantly, these over- 
laps were highly statistically significant whereas genes that 
were reciprocally regulated upon BAPl and FOXK2 de- 
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Figure 5. FOXK2 and BAPl regulate a common set of target genes. (A) Venn diagrams showing the overlap genes which are commonly down-regulated 
(top) or up-regulated (bottom) following depletion of FOXK2 (12) or BAPl (22) in U20S cells. (B) Comparison of the numbers of overlapping genes 
showing up- or down-regulation following depletion of FOXK2 or BAPl (blue bars) compared to the numbers expected by chance (grey bars). P-values 
are shown above the columns. (C) Validation of FOXK2 and BAPl co-regulated genes. FOXK2 or BAPl was depleted in U20S cells and the indicated 
target gene expression was detected by RT-qPCR. Non-targeting (NT) siRNAs were used as a control. Data are shown relative to the expression seen with 
NT siRNA (taken as 1) and are the averages plus standard deviations (error bars) from three independent experiments. 



pletion were far fewer in number and were either insignif- 
icant or lowly significant compared to randomized datasets 
(Figure 5B). This indicates a common mode of function 
for FOXK2 and BAPl in either up- or down-regulating 
gene transcription. To focus on potential direct FOXK2 tar- 
gets we associated the commonly deregulated genes with 
FOXK2 binding events (within ±5 kb from the TSS) from 
our ChlP-seq data, leaving 184 (74%) commonly down- 
regulated and 263 (71%) commonly up-regulated genes. We 
selected one gene from each class of co-regulated genes 
for detailed verification and subsequent further analysis; 
TP 5 3 13 (up-regulated upon depletion) and H2AFX (down- 
regulated upon depletion). Depletion of either FOXK2 or 
BAPl led to enhanced TP 5 3 13 expression whereas H2AFX 
expression was decreased as expected from the microarray 
analysis (Figure 5C). We also tested combinatorial deple- 
tion of FOXK2 and BAPl, and while increased activation 
of TP53I3 was observed, no further reductions in H2AFX 
occurred upon depletion of both FOXK2 and BAPl (Sup- 
plementary Figure S5). 

Next, we verified FOXK2 and BAPl occupancy in the 
regulatory regions of these genes, and binding of both fac- 
tors was detected (Figure 6A). Upon FOXK2 depletion (e.g. 
Supplementary Figure S6B), both FOXK2 and BAPl bind- 
ing was diminished at both the TP53I3 and H2AFX loci, 
indicating a role for FOXK2 in BAPl recruitment (Figure 
6A). To further extend this finding, we also tested two re- 
gions that were identified by our ChlP-seq analysis of Flag- 
tagged BAPl. We verified BAPl binding and both regions 
also exhibited FOXK2 occupancy. Importantly, FOXK2 
depletion also caused reduced BAPl binding to these loci, 
further emphasizing its role in BAPl recruitment (Figure 
6 A, bottom graphs). We also tested whether there was a re- 



ciprocal effect on FOXK2 binding upon BAPl depletion 
but little effect was seen on the ChIP signal for FOXK2 
following reductions in BAPl expression levels by siRNA 
treatment (Figure 6B; Supplementary Figure S6B). 

BAPl has been shown to be responsible for his- 
tone H2AK119 deubiquitination (32), therefore we tested 
whether BAPl functioned in this manner at the regulatory 
loci controlled by FOXK2 and BAPl. Depletion of BAPl 
caused increases of histone ubiquitination at all of the reg- 
ulatory loci, albeit often to moderate levels (Figure 6C). 
Depletion of FOXK2 had less effect, but depletion of both 
factors together caused a significant increase in H2AK119 
ubiquitination at all loci tested (Figure 6C). Longer term 
depletion of FOXK2 in stable FOXK2-depleted cell lines 
led to bigger increases in H2AK119 ubiquitination at the 
TP 5 3 13 and H2AFX regulatory regions, consistent with its 
role in BAPl recruitment to these loci (Supplementary Fig- 
ure S7). 

Collectively, these results demonstrate that FOXK2 func- 
tions together with BAPl to control transcription in either 
a positive or negative manner. In this context, FOXK2 acts 
to recruit the PR-DUB complex component BAPl to chro- 
matin and thereby promote local histone deubiquitination. 

DISCUSSION 

Transcriptional control of gene expression is controlled in 
part through localized changes in chromatin architecture. 
This is driven by a wide array of different histone modi- 
fications which themselves are generated under the direc- 
tion of sequence-specific transcription factors and are of- 
ten dynamically modified in response to changes in the cel- 
lular environment (reviewed in (33)). Here, we have stud- 
ied the forkhead transcription factor FOXK2 and demon- 
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Figure 6. FOXK2 acts as a chromatin targeting factor for the BAP 1 -containing deubiquitinase complex. ChIP analysis of endogenous FOXK2 and BAPl 
binding to the indicated genomic loci in U20S cells treated with non-targeting control siRNAs (siNT, black bars) or siRNA against FOXK2 (siFOXK2, 
white bars) (A) or siRNA against BAPl (siBAPI, grey bars) (B). ChIP was performed with nonspecific IgG or anti-FOXK2 and BAPl antibodies; the 
data are shown relative to the input DNA and are the averages of two (B) or three (A) experiments. (C) ChIP analysis of H2A and H2AK1 19 ubiquitin 
levels at the indicated genomic loci in U20S cells treated with non-targetting control siRNAs (siNT; Ught grey bars), siRNA against BAPl (siBAPI; dark 
grey bars), siRNA against FOXK2 (siFOXK2; white bars) or both BAPl and FOXK2 (black bars). ChIP was performed with non-specific IgG, anti-H2A 
and ubiquitinated K119 H2A antibodies (H2Aub); the data are the averages of three experiments. Statistically significant differences between siNT and 
siBAPI /siFOXK2 treated samples are shown; P-value <0.05 (*) and <0.01(**). (D) Model depicting the role of FOXK2 in nucleating the recruitment 
of chromatin remodelling complexes to chromatin. FOXK2 recruits both the PR-DUB and the SIN3A complex, potentially through direct or indirect 
interactions with the shared subunit HCFCl. Question marks denote uncertainty about which interactions are direct. BAPl can then cause local histone 
deubiquitination, and histone deacetylation can potentially be achieved through the HDACl component of the SIN3A complex. 



strate that it can recruit the PR-DUB complex to chro- 
matin, where it can modify the local chromatin environment 
by promoting histone deubiquitination through its BAPl 
subunit (see model in Figure 6D). While we have provided 
most evidence to support this element of our model, we also 
demonstrate that FOXK2 can also bind to the SIN3A co- 
repressor complex, potentially allowing localized histone 
deacetylation as well as deubiquitination. 

BAPl is becoming increasingly recognized as an impor- 
tant protein in the context of tumourigenesis. BAPl is often 
mutated or deleted in a range of cancers, most commonly 
in uveal melanomas, mahgnant pleural mesothehomas and 
clear cell renal carcinomas (reviewed in (29,30)) indicating 
that it represents a tumour suppressor protein. Recently, 
its importance was further emphasized as an extensive can- 



cer sequencing study found BAPl to be one of seven genes 
whose mutation is the most indicative of poor prognostic 
outcome (30). Previous studies have suggested that HCFCl 
binds to BAPl and can therefore potentially recruit it to 
chromatin (28,34). As HCFCl can also bind to E2F tran- 
scription factors, this provides a potential route to specific 
recruitment of BAPl to regulatory regions in chromatin 
(35). Indeed, BAPl depletion leads to changes in E2F tar- 
get gene expression, although a role for E2F in nucleating 
this event was not demonstrated (36). More recently, YYl 
was shown to permit recruitment of BAPl to the regulatory 
region of the C0X7C gene to activate its expression (22). 
Here, we provide an alternative route through which BAPl 
can be recruited to specific genomic regions through the 
transcription factor FOXK2. As FOXK2 occupies thou- 
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sands of genomic loci, many of which are found in open 
chromatin (12), this provides a molecular mechanism to 
help explain how BAPl might function across the genome. 

In addition to binding to the PR-DUB complex, we also 
found that FOXK2 can bind to the SIN3A co-repressor 
complex, with interactions detected with SINS A, RBBP4, 
RBBP7 and HDACl. This is in keeping with the obser- 
vation that the closely related FOXKl can also bind to 
the SINS complex (7,8), and suggests a repressive role for 
FOXK2 in this context. One way in which FOXK2 might 
repress transcription is through the HDACl subunit of this 
complex, which would lead to localized histone deacetyla- 
tion and hence changes in chromatin structure. It is possible 
that this activity might synergize with the effects of histone 
deubiquitination through the PR-DUB complex, but it is 
unclear whether FOXK2 recruits these to the same regula- 
tory loci either at the same time or in a sequential manner. 
However, it is intriguing that the FOXK2 interacting pro- 
tein HCFCl has previously been shown to also bind to both 
SIN3A and BAPl (Figure IC; (28,34,37)), thereby provid- 
ing a potential platform for recruiting both complexes to 
FOXK2 (see model in Figure 6D). However, it is unclear 
whether FOXK2 makes direct interactions with HCFCl, 
BAPl, subunits of the SIN3A complex or combinations of 
these. 

Histone H2AK1 19 ubiquitination is typically thought of 
as repressive in nature (38), and thus the recruitment of 
the BAPl deubiquitinating enzyme would be expected to 
promote transcriptional activation. However, the PR-DUB 
complex is not generally associated with transcriptional ac- 
tivation mechanisms and instead appears to be repressive 
in nature, potentially through promoting histone ubiquiti- 
nation dynamics (reviewed in (39)). Indeed, this is consis- 
tent with our demonstration that BAPl depletion causes 
increased ubiquitination levels at regulatory regions associ- 
ated with genes that are either repressed or activated under 
the same conditions (see Figures 5C and 6C). It remains un- 
clear though whether the changes in histone ubiquitination 
we observe have a direct influence on the outcome of tran- 
scription at associated target gene loci. For example, tran- 
sient depletion of FOXK2 has little effects on steady state 
histone ubiquitination levels (Figure 6C) and yet influences 
both TP53I3 and ^2 v4FX transcription (Figure 5C). Fur- 
thermore, FOXK2 and BAPl depletion synergizes to cause 
enhanced TP53I3 gene activation (Supplementary Figure 
S5), but no synergy is seen at the H2A ubiquitination level 
(Figure 6C). Thus, although it is clear that FOXK2 recruits 
BAPl and influences histone ubiquitination levels, there is 
no clear link between the chromatin changes and gene reg- 
ulation. These observations point to additional important 
regulatory activities of FOXK2 and/ or BAPl which may 
function in a locus-specific manner. Moreover, it is possible 
that FOXK2 might coordinate the recruitment of repres- 
sive and activating complexes to the same loci to turn the 
regulatory regions on or off under different conditions. Al- 
ternatively, FOXK2 might be repressive in nature at some 
loci but activating at others. Indeed, in keeping with this 
hypothesis, we previously demonstrated that roughly equiv- 
alent numbers of direct FOXK2 target genes were up- and 
down-regulated following FOXK2 depletion (12). 



Our demonstration that FOXK2 can function in modify- 
ing chromatin structure through recruiting histone deacety- 
lation and deubiquitination complexes, adds to a growing 
number of connections with chromatin-associated events. 
For example, one recent study identified FOXK2 as a com- 
ponent of 5-hydroxymethylcytosine-DNA binding protein 
complexes (40) whereas another studied indicated that 
FOXK2 showed a preference for binding 5-formylcytosine 
(41), suggesting that there is a relationship between FOXK2 
and regions containing either methylated DNA or its 
derivatives. More recently, we have shown that FOXK2 and 
BAPl are also components of complexes containing the 
MBD5/6 methyl binding domain family members (42). It 
will be interesting to determine whether FOXK2-mediated 
recruitment of the PR-DUB complex has functionally im- 
portant interactions with other chromatin components that 
recognize the underlying methylation state of the genome. 

SUPPLEMENTARY DATA 

Supplementary Data is available at NAR online. 
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